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We study the evolution of the cosmological parameters, namely, the deceleration parameter q{z) 
and the parameter of effective equation of state in a universe containing, besides the ordinary matter 
and dark energy, a self-interacting (collisional matter), in the generalized f{R,T) theory of gravity, 
where R and T are the curvature scalar and the trace of the energy-momentum tensor, respectively. 

We use the generalized FRW equations, the equation of continuity and obtain a differential equation 
in H{z), and solve it numerically for studying the evolution of the cosmological parameters. Two 
f{R,T) models are considered and the results with collisional matter are compared with the ones 
of the ACDM model, and also with the model where there exists only non-collisional matter. The 
curves show that the models are acceptable because the values found for Wsff are consistent with 
the observational data. 

PACS numbers: 04.50.Kd; 98.80.-k; 95.36.-l-x 


I. INTRODUCTION 

It is well known nowadays that our current universe is experiencing an accelerated expansion [l]-[7]. There are 
two general ways to explain this accelerated expansion of the universe [8]-[12]. The first way is considering that the 
universe is essentially filled by an exotic fluid with negative pressure, responsible of it acceleration, called the dark 
energy. The second way is modifying the gravitational action from the General Relativity (GR) without the need of 
the dark energy, firstly considered by 1970 [13]. Buchdahl has generalized the Einstein equations by substituting the 
Ricci’s scalar R by an arbitrary function of R, named /(i?) theory of gravity. Another theories, view as alternative 
to the GR, also have been undertaken, still in the way to better explain the obscure content of the universe, as, 
/(i?, G)[14]-[15], /(i?, T)[16]-[17] and /(i?,T, i?^yT'^^)[18]-[19], where G, T, R^^ and are the invariant of Gauss- 
Bonnet, the trace of energy-momentum tensor, the Ricci’s tensor and the energy-momentum tensor corresponding to 
the ordinary content of the universe. 

In this paper, we focus our attention on f{R,T) theory of gravity. This theory has been considered first by Harko 
and collaborators [20]. Another authors also have considered this present theory and important results have been 
found [21]- [29]. 

However, any one of the works developed in these papers does not solve the coincidence problem, that is, how does 
the universe transits from the decelerated phase to the current accelerated phase? [31]. The universe is accelerated 
for the low values of the redshift and the one that splits the two phases (transition redshift), denoted Zt, those the 
current value is Zt = 0.46 ± 0.13 [32], [33]. The f{R,T) theory can successfully predict the transition from the 
matter dominated phase to the current accelerated one and several works have described this transition [34]. The 
reconstruction of f{R, T) models describing the matter dominated and accelerated phases had been performed in [30]. 
Various works also have been performed, still in the optic to better explore this transition with interesting results( 
see [35], [36] and [43]). 

In this paper, we focus our attention on the epoch after the recombination where, beside the well known ordinary 
matter (the dust), there is the self-interacting matter, called collisional matter. The model of collisional matter has 
been studied in some works within others theories of gravity, leading to interesting results [38],[44], [45], [46]. This 
approach of considering new form of matter, instead of the cold dark matter can edify us on the choice of the models 
of modified gravity. Oikonomou and collaborators [40] have studied the time evolution of the cosmological parameters 
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during the late time acceleration of the universe with the presence of the collisional matter in the framework of modified 
f{R) gravity. In this paper, we have extended the same idea to the f{R, T). Some /(!?, T) models have been considered 
and the behaviours of the cosmological parameters have been performed and compared with the ACDM model. We 
see within many results that depending on the input parameters according to the model under consideration, the 
inclusion of collisional matter may lead to a better explanation of the phase transition, comparatively to the model 
where just the usual ordinary matter is considered. 

The paper is organized as follows: In section II we describe the general formalism of /(i?, T) theory of gravity. The 
collisional matter that self-interacts is presented in the section III. The section IV is devoted to the study of time 
evolution of the cosmological parameters where the universe is considered to be filled by the usual ordinary matter 
and the collisional one. Here, for the study of these cosmological parameters, we focus our attention on the transition 
from the decelerated phase to the accelerated one. In the section V we examine the evolution of the equation of 
state of the dark energy where the matter content is assumed as a fluid is composed by the collisional matter and the 
radiation. The conclusion and perspectives are presented in the section VI. 


II. GENERAL FORMALISM IN /(R, T) GRAVITY 


In this section we present the generality of f{R,T) theory by substituting the curvature scalar R of the GR by a 
function of R and the trace T, and writing the action as [37] 


S= 




( 1 ) 


where i?, T denote the curvature scalar and the trace of the energy-momentum tensor, respectively, and = SttQ, Q 
being the gravitation constant. 

The energy-momentum tensor associated to the matter is defined by 


T — ^ d {y/ gCm) 

^ ^ 

Let us assume that the matter Lagrangian density Lm only depends on the components of the metric tensor, but not 
on its derivatives. Thereby, one gets 


T 

fliy 


g^i'R’n 


2dL^ 


( 3 ) 


Within the metric formalism, varying the action (I) with respect to the metric, one obtains the following field equations 
[37] 


IrRiiu - -^giJ.i^f{R,T) + u)fR = - friTyi, + Qy^) , 


( 4 ) 


where fa, fr are the partial derivatives of f{R, T) with respect to Rand T respectively. The tensor Qyi, is determined 
by 


Qy. = 9° 


ajS STgp _ 




— —2Tyi, -|- gyi/Cm — 2g' 


a/3_ 




dg>^''dg°‘d 


( 5 ) 


As mentioned in our introduction, we assume that the whole content of the universe is a perfect fluid. Then, by 
setting the matter Lagrangian density to Cm = —Pm, the energy-momentum tensor may take the following expression 
(see appendix) 


Cyi/ — (Pm T Pm) RyRu PmPyi^^ (6) 

where the four-velocity satisfies the relations UyU^ = I and u^V^Uy =0. Note that the expression (6) is obtained 
with the consideration according to what the pressure does not depend on the metric tensor. Within 
this consideration according to what the Lagrangian density does not depend of the metric tensor, the 
contribution of the last term of (5) vanishes and this eqnation takes the following form (see also Eq. 
(21) of [20] ) 


^ yy — 2Tyi/ PmPyu- 


( 7 ) 
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However, in the case of the electromagnetic field or scalar field, there is actively contribution of the 
last term of (5) because in these cases the respective Lagrangian densities £ = g°‘"'gPcrFapF^cr and 
C = gf^‘'dfj_4’diy(l), explicitly depend on the metric. (See also comments after Eq.(18) and Eq.(19) of [20]). 

Making use of Eq (7), the field equation (4) can be reformulated as 

~ ^ (8) 
where Kgjy = is the effective gravitational constant and 


T^e// 

fliy 


T, 




/t (2 


gtiuif RJr) F fTPingfiu 


(9) 


representing the effective energy-momentum tensor of the matter. 

Taking into account the covariant divergence of Eq.(8), with we use the identities — \Rgfj,v) = 0, and 

Jr = -Rf,i,Vf^fR, on gets 


V^T —_ 

V -L 111/ - 2 J- 

- It 




fiiy 


[Til, + Qii.) V'^/t - ^gnufrV^T 


After some elementary transformations, Eq. (10) takes the following expression 


( 10 ) 


= 


It 


- ^^Kfrdiip + (1 - Sw) 




( 11 ) 


It comes from the above expression that for the /(i?, T) theory of gravity, the conservation of the energy-momentum 
tensor is not automatic. In order the reach this condition, it is important to constraint the second hand of this 
equation to zero (see Eq. (18)). This condition is extremely crucial for determining the suitable expression of the 
Trace depending part of the algebraic action. The reference [51] considers the constraint of the conservation 
to energy-momentum tensor (as in (18)) and one obtains, after the integration on the right side of 
(11), /(T) = + Pot with ao,f3o real constants; (see (19) of [51]). 

Here we are interested to the spatially flat Eriedmann-Robertson-Walker (FRW) metric 


ds'^ = dt'^ — a{t)'^[dx'^ -I- dy'^ + dz'^]. 


( 12 ) 


Thus, the Ricci scalar in this background is given by 

R= + H). 


(13) 


III. COLLISIONAL MATTER MODEL WITHIN f(R, T) THEORY 

This king of matter has been introduced firstly in [38, 39] in the framework of GR and f{R). The evolution of the 
collisional matter depends on the source that drives the whole content (matter and energy) of the univers. 

The study of the model of collisional matter with perfect fluid has been performed widely in [40]. In this work, we 
recall that our basic assumption is that the matter has a total mass-energy density, denoted by £m, which is assumed 
to be depending on two contributions as 

— Pm F PmR- ( 1 ^) 

Here, pm refers to the part that does not change due to the fact of characterising the usual matter content [38], [39]. 
Concerning the term PmR, it expresses the energy density part of the energy momentum tensor associated with 
thermodynamical content of the collisional matter. This fluid could not be dust, but possesses a positive pressure and 
satisfying the following equation of state 

Pm — R^Pm^ ( 1 ^) 

where w denotes the parameter of equation of state of the collisional matter whose values is such that 0 < w < 1 . 
The potential energy density is assumed to have the form [38],[39] , 

n = Ho Fwln{-^). 

PmO 


( 16 ) 
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where the constants pmo and TTq are their current values. 

Accordingly, after elementary transformation the total-energy density of the universe can be written as 

Em = Pm(^ + no + wln{-^)\ (17) 

\ PmO ) 

Through the continuity equation, the motions of the volume elements in the interior of a continuous medium can be 
traduced 


= 0, (18) 
and the energy-momentum tensor takes the following form 

where = dx^/ds is the four velocity, satisfying the relation u^Uv = 1. We note here that Pm = Pm because 
the pressure of ordinary matter is negligible. Making use of the FRW line element (12), the conservation law of the 
Equation (18) 


Em “b 3 \Em “b Pm) — O 5 
a 

which, upon consideration of the equations (15) and (17), leads to [38] 



( 20 ) 


( 21 ) 


where oq is the current scale factor. 

We can describe the collisional matter by equations (17) and (21) and use therm in the rest of the manuscript in 
the next section. The value of Hq is equal to [38] 

We will use the same value [38],[42] of TTg = 2.58423 for making the numerical study. According to the fact that 
p{to) = po, and making use of Eq. (16), one gets 77(to) = no, meaning that this latter is the current value of the 
potential energy per unit rest-mass, associated with the infinitesimal deformations (expansions or/and compressions) 
of the fluid (see detail in [39]). 


IV. LATE TIME COSMOLOGICAL EVOLUTION IN f{R,T) THEORY 
A. Deceleration Parameter 


In this subsection we examine the deceleration parameter q{z) in f(R,T) considering a that the universe, beside 
the ordinary matter and dark energy, also is filled with collisional matter at late time cosmological evolution. In this 
optic, we rewritte the field equation (4) in the following form 


3iJ2 


1 + /t 



1 

Ir 


l{f-RfR)-^RHfnR+p^fT 


(23) 


-2H - 377^ 


1 + /r 
Ir 


1 

Pm “b ~ 
JR 


2777?/flfi + R/rr + R^/rrr 


- RIr) - Pm fr 


(24) 


where we have assumed = 1. The parameter H = ^ denotes the Hubble parameter, the “dot”, the time derivative 
of d/dt- We propose to reformulate the right side of (23) and (24) respectively in term of effective energy density peff 
and pressure Peff, and we write (23) and (24) as 




3H^ = rI 


( 25 ) 
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Here, the effective energy density peff and effective pressure Peff are respectively defined as 


P^ff ~ 


1 + /t 


{f-RfR)-iRHfRR+PmfT 


(26) 


(27) 


Peff — Pm “t“ 


1 

1 + /t 


‘2‘HRfRR + RJrr + R^Irrr 


- RIr) - Pmfr 


(28) 


In this way, we have a matter fluid representation of the so-called geometrical dark energy in f{R,T) gravity with 
the energy density poE = Peff - £m and pressure poE = Peff - Pm ■ 

From the conservation law, the effective energy density evolves as 

By making use of the equations (23) and (24), (29) can be rewritten as 

fRB •• • ■ n 1 -I- fr fr f 

+ 4HH) + 3{H + H^) + ^^Sm+Pm^ + T^ = ^- (30) 

JR JR JR 2Jr 

One can rewrite the above equation using the redshift z = i — 1 as follows 


cPH 3 dH I UhV 3/fi(^g^-(l + z)H^^ +|+p^/r + (l + /r)£m 
dz2 1 + z dz H\dz ) ld,H^fRR(l + zf 


(31) 


As convention in this paper, we set to unity the current value of the scale factor. We also assume the algebraic function 
as a sum of two independent functions, such that, /(i?,T) = f{R) + f(T). The f(T) model is the one obtained by 
imposing the conservation of the energy-momentum tensor (11). 

The dynamics and stability of this model are studied in [37] and interesting results have been found. This model 
reads 


/(T) cx r“, 

with a = 2 {i+w) ■ Here the trace of the energy- momentum tensor (19) reads T = Sm — ^Pm- 
From Equations (17) and (21), the model f(T) and its first derivative are 


/(T) = (1 + z)3“p„o 


1 -I- ilo — 3w(l — ln(l + z)) 


(32) 


(33) 


/T=a(l+z)'^““'VmO^“"'^ 


1 -I- IIq — 3r(;(l — ln(l -b z)) 


(34) 


By using (33) and (34), one rewrite (31) as 


d^H 3 dH 1 fdHV 

- - Ji(z) - J2(z) - J3(z), 


1 -b z dz H \ dz J 


(35) 


where 


Mz) = 


ISH^fRRil + zf 


(36) 


J2{z) = 


Pmo(l + 2 ) ( 1 + Ro + 3rc ln(l -b z) 
l^H^fRR 


(37) 
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J3(z) = ^a(l + w) + i(l - ^(1 + 2;)^“ ^ ^1 + iTo - 3u;(l - lii(l + z))^ J2(z) (38) 

The equation (35) is the main tool to be used for analysing the cosmological evolution of the considered f{R,T) 
models. The numerical analysis also will be performed in order to find the Hubble parameter. In a first step, by using 
the expression 


q{z) 


(1 + z) dH 
H{z) dz 


(39) 


we can perform the study of the deceleration parameter with the account of the collisional matter, this latter assumed 
as a function of the redshift z. In this rubric we will focus our attention on the transition from the decelerated 
phase to the accelerated one. Therefore, we will compare the results with the case where the non-collisional matter 
(pressure-less) is considered, and also compare with the ACDM case. 

We recall that in [42], the corresponding deceleration parameter of the expansion in the flat friedmann cosmology 
for the ACDM model is given by 


q{z) 


~Y~ 


(1 + 2 :) — Lq 


f^mo(l + 2)^ + Lq 


-1 


(40) 


where flmo = 0.279 and Lq = 0.721 [42]. 

Moreover, the effective equation of state parameter Wg// is defined as, 

(41) 

(42) 

and we present its evolution versus the redshift. A comparison with the results to be found will be done with the 
results coming from pure f{R,T) models with non-collisional matter. 


Weff = 


Peff 

Peff 


= -1-b 


2(1 -1- z) dH 
3H{z) dz ’ 


B. Study of particular cases f{R) models in modified f{R,T) gravity universe with collisional matter 

In this subsection we discuss two f{R) models, performed in [40, 43] for describing the late time cosmological 
evolution of our universe under the consideration that it is filled with collisional (self-interacting) matter. Thereby, 
we solve numerically the differential equation about H(z). 

The task here is to find the models that should be in agreement with the observational data, namely the deceleration 
parameter q{z) and the parameter of effective equation of state tCg//. 

The first model considered is a modified power-law model of the curvature scalar 

/(i?)=Ao(A + i?)"; (43) 

where Aq, A and n are positive constants. 

The second, is an exponential model 

f{R) = Roe^^, (44) 

where Rq and P are constant parameters. 

Numerical plots of the deceleration parameter g(z) and the effective equation of state Weff are presented from (35) 
with collisional matter {w = 0.6) and non-collisional matter (w = 0) for each model. 

For Fig. 1, related to the first f{R,T) model, we see that, at the left hand side, traducing the evolution of the 
parameter of deceleration the curve representative of the collisional matter is nearby that of ACDM than the curve 
representative of the non-collisional matter. We also see that, in the order, collisional matter, ACDM and non- 
collisional matter, the transition from the decelerated phase to the accelerated phase is realized from high to the low 
value of the redshift. Concerning the curve representative of Weff the same behaviours of the phase transition order is 
confirmed. We see that the behavior of the cosmological evolution of the deceleration parameter q{z) and the effective 
equation of state parameter Weff for the model f{R,T) = Ao(A -|- i?) -|- with and without collisional matter are 
similar to the model f{R) = Ao(A-|-i?)" without collisional matter obtained by the authors of ([40]). Moreover, as 
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Figure 1: The graphs show the evolution of q{z) and Weff versus z for the model f{R, T) = Ao(A + Rf" + T“. The 
graphs are plotted for Aq = 1, A = 13.5 and n = 0.5 The red, blue and the magenta refer to non-collisional, 

collisional matter and ACDM model, respectively. 




Figure 2: The graph shows the evolution of q{z) and Weff versus z for the model f{R, T) = R^e^^ + T°‘. The graph 
are plotted for i?o = 1 and /3 = 1.5. The red, blue and the magenta refer to non-collisional matter, collisional matter 

and ACDM model, respectively. 


the same authors, we see that the transition from the decelerated phase to the accelerated phase is realized from 
high to the low value of the redshift in same order of the content of universe. Also, we remark that, the transition 
of the redshift obtained by this authors for f{R) model without collisional matter is equal to the one obtained for 
f{R,T) model without collisional matter, whereas the transition of the redshift to the f{R) model with collisional 
matter is strongly than the transition of the redshift for f{R,T) model. For this observation, we can conclude that 
the depending of the trace T of the energy-momentum tensor in this model is very important to reproduce the ACDM 
model. 

Regarding the second model, whose cosmological parameters are plotted at Fig. 2, we see that for the transition from 
the decelerated to the accelerated phases, from the high to low redshifts, the order is non-collisional matter, ACDM 
and collisional matter. Moreover, we see that in the accelerated phase, the collisional matter is more nearby the 
ACDM than the non-collisional one, while in the decelerated phase, both the collisional and non-collisional matters, 
almost follow the same trajectory for high redshifts. This aspect is very important because revels that the collisional 












takes its origin from the ordinary matter. We then conclude that at early times there is no collisional matter, but 
just non-collisional, and as the time evolves, the collisional matter start being created from the non-collisional matter. 
After crossing the transition line, the collisional matter approaches more the ACDM as the low redshifts are being 
reached. We see that the behaviour of the cosmological evolution of the deceleration parameter q{z) and the effective 
equation of state parameter Weff for the model f{R, T) = + T“ with and without collisional matter is similar 

to the model ,f{R) = Rqc^^ with and without collisional matter obtained by the authors of ([40]). At the high value 
of the redshifts, the trajectory of this parameters is bewildered, both the collisional and non-collisional matter for 
f{R) and f{R,T) models 


C. f{R,T) Models with Cardassian Matter 

The cardassian model is firstly developed by [44],[45],[46] and is another approach to describe the self interacting 
matter. According to the observational it is well known that the universe is flat and accelerating; the cardassian 
model of matter is characterised by negative pressure and there is no vacuum energy whatsoever. In this subsection 
we perform the same study as in the previous section and use the models f{R) models (43), (44). According to the 
so-called Cardassian model of matter [44] , [45] , [46] , the total energy density Em of matter is defined by 

em = P + pK{p) (45) 

where p is related to the ordinary matter-energy density and pK [p) describes the cardassian model of matter and is 
general, a function of the ordinary mass-energy density p. 

According to the original Cardassian model, one has 

pK{p) = Bp^\ (46) 

where B is a real number and it is fixed n' < | in order to guarantee the acceleration. 

By taking into account the relation (46), we reobtain the total energy density of matter (45) as, 

e^ = p + Bp^'. (47) 

It is important to point out here that the above expression may be expressed in term of the redshift z. To do so we 

take into account (21) and cast it on the form p = Pmo(X + ^ where we assumed the current scalar factor oq = 1. 

Therefore, the equation (47) becomes 


Em — Pmo(l + ^)^ ^^(l-|-z)^^ 


(48) 


Following the same idea as [40], we assume that the late time evolution is governed by the geometric dark fluid 
with negative pressure coming from the f{R,T) models plus a gravitating fluid of positive pressure, satisfying the 
following equation of state 


P = Wkp. 


(49) 


We mention in this model that the pressure has to be assumed as negative quantity, p < 0, because the cardassian 
model of matter possesses negative pressure. In such a situation, this negative has to be assumed as the one which 
generates the acceleration of our universe. 

Adopting the assumption that the total energy density of matter is (47), the main equation (31) can be used for 
analysing the evolution of the viable f{R,T) models considered with Cardassian self-interacting matter, yielding 


(fH 

~d^ 


3 dH 
1 + z dz 


-J[{z)-J'^{z)-J'^{zl 


(50) 


where 


18ff3/^fl(l + zf 


(51) 
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J2{Z) 


Pmo(l + ^) + 


(52) 


) a —1 

J'(z). (53) 

By adopting the same treatment as in the previous section, we present the evolutions of the deceleration parameter 
q{z) and the effective equation We// of state for the late time cosmological evolution of the universe with the /(i?, T) 
models including Cardassian Matter. 


— PmO°‘ ^(1 + z)^^ \a + 1/2)(1 + Wk)^ + Wk + ^^(1 + z)^^ ^ 




Figure 3: The graph shows the evolution of q{z) and Wgff versus z for the model f{R, T) = Ao(A + i?)" + T“. The 
graph are plotted for Aq = 1, A = 13.5 and n = 0.5 The red color, blue and the magenta refer to collision-less, 
cardassian and ACDM model, respectively and we use B = 0.2, n' = —3. 

Here, for the first model. Fig. 3, we see that the transition from the phase of deceleration to the acceleration one 
is realized, from the high to low-redshifts, in the order, ACDM, cardassian matter and collision-less matter. In the 
decelerated phase, both the cardassian and collision-less matters are confused, while, within the accelerated phase, the 
cardassian matter detached from the collision-less one and start approaching the ACDM model as the low-redshifts 
are being reached. As the authors of ([40]), we see, both the cardassian matter and collision-less matter, that the 
cosmological evolution of the deceleration parameter q(z) and the effective equation of state parameter Weff for the 
models f{R,T) = Ao(A-|-i?)" -|- T°‘ and /(i?) = Ao(A-|-i?)" are similar behavior and are same trajectory at high 
value of the redshifts. The transition of the redshift for each model with cardassian matter and collision-less matter 
is inferior at the transition of the ACDM model. 

Concerning the second model. Fig. 4, the cardassian and collision-less matter are also confused for high redshifts 
and as the low-redshifts are being reached, the cardassian matter approximate the ACDM model. We observe at the 
low value of the redshifts, with cardassian matter, that the behavior of the f[R^ T) = Roe^^ + T°' model is not similar 
of the model f{R) = Rqb^^ obtained by the authors of ([40]), but, when the redshift 2 increase, the trajectory of this 
models with cardassian matter are similar behaviour and bewildered with the collision-less matter trajectory. 


V. EQUATION OF STATE FOR DARK ENERGY WITH COLLISIONAL MATTER FLUID IN f{R, T) 

GRAVITY 

The work performed in this section is about the cosmological evolution of the universe through f{R,T) models in 
presence of matter composed by collisional matter and radiation. 
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Figure 4: The graphs show the evolution of q{z) and Weff versus z for the model f{R,T) = + r“. The 

graphs are plotted for Rq = 0.5 and /3 = 1.4. The red, blue and magenta refer to collision-less matter, cardassian 
matter and ACDM model, respectively and we use B = 0.2, n' = —3. 


We focus our attention on the exponential f{R) model [47],[48] given by 

f{R) = R-PRsil-e-^/^n, ( 54 ) 

where ci = —f5Rs and C 2 = Rs constants with curvature scalar dimension. The conditions for having a viable 
exponential model has been studied in [49]. 

Within flat FRW metric and considering the above content of the universe, the field equations, from the general 
form (4), take the following form 

= (1 + fT)Pmatt + \if- R/r) - ^HJr + Pmatth. (55) 

- 2fRH = (1 + fT){Pmatt + Pmatt) + /fl " (56) 

where pmatt and Pmatt denote the energy density and pressure of all perfect fluids of generic matter, respectively. 

Concerning the component to be considered here, we assume that the universe is filled with collisional matter 
(self-interacting matter) and the relativistic matter (radiation). This means that in our analysis, the contribution 
from the radiation may play an important role. Therefore, the matter energy density, Pmatt, is given by 

Pmatt — PrOR : (57) 

where pro is the current energy density of radiation. The pressure of all perfect fluids of matter is given by 

Pmatt ~ Pm T Pr- (53) 

The first term characterises the pressure of collisional matter, given by (15), and the second, is related to the radiation. 
Making use of (17) and (21), one can put Eq. (57) into the following form 

Pmatt = Pmoa~^ i 1 + Uq + 3w ln(a) ) -f Proa~‘*^. (59) 


In the same way, Eq. (55) can be rewritten as 


H2 _|_ -(X _ -1_ + ^(p„ 

alnab dlna 3 3 


4 “ Pmatt ); 


(60) 
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while the scalar curvature R can be expressed as 


dH 

R=-( 12 i ^2 + 6 i^—-). (61) 

a In a 

Through the Eq. (59), we assume that the total matter-energy density becomes 

Praatt = Pmoig{a) + Xa"'‘)- (62) 

In the previous expression, y is defined by y = pro/pmo — 3.1 x 10“'^, pro being the current energy density of radiation. 
The parameter g{a), describing the nature of the collisional matter (view as perfect fluid), is equal to 


g{a) = a I 1 -I- TTq -I- 3w ln(a) I. 


(63) 


Note that for the non-collisional matter (assumed as the dust), for which the parameter w = 0, one gets g{a) = a~^. 

In the optic to better study the cosmological evolution of the f{R,T) models in the framework of fiat FLRW 
universe, we may introduce the variable [41] 


Pde , , _4 

VH = -= — - g[a) - ya , 

Pmo 

R dq(a) 

VR = — - -7\ -> 

d In a 


(64) 

(65) 


K PmO 


where poE denotes the energy density of dark energy, rrr, the mass scale, given by [50] mr = g 
Making use of (60), the expression ^^ 37 ^ yields 


1 


1 dR 


m? din a fun 


fr 1 1 

{Pmatt 4“ Rmatt^ 4“ _ n Pmatt _ n 4“ ^ _ 9 (/ R) (1 f r) 

6m^ m Dm 


3m^ 


H dH id" 


m? din a m^ 


■ ( 66 ) 


We recall in this work that the trace of the stress tensor depends on the nature of the matter content. Therefore, 
using the same model (32) where T = pmo[fl(a)(l ~ 3r(;)], and also (58) and (59), Eq. ( 66 ) becomes 


1 dR 


m? din a H^/rr 


3^m2^(l -b /3) ^g(a)(l - 3w)^ ^(1 -b w)g{a) + ^ya - yn + 



dH 


\m'^ 

din a 

m^J 


where f3 = a — 1. Eqs. (61) and (67) are reduced to a coupled set of ordinary differential 


dyH 1 , 4dg(a) ^ , 

-7] -= -oVR - - q Ti- 

dm a o 6 dma 


(67) 


( 68 ) 


dyR 

din a 


<Pg{a) 

d In 


3'^m2^(l-b^)^5(a)(l - 3«;)^ ^(1 4-w)5(a)-b ^ya 
, 1 .. pN . n/^1 dyu , ld5f(a) 


iRRm^iyH + g{a) + XO' '^) 


Moreover, the curvature scalar is expressed as 


R = —3m^ ( dj/ij -b 45 (a) 4- 


dyn , dg{a) 


d In a d In a 


(69) 


(70) 


By operating the differentiation of the relation ( 68 ) with respect to In a, and eliminating from this result and 
(69), we obtain 


d‘^y 


H 


din a 


+ 4- 


fR 


dyn 


fR-2 


yn 


d^g{a) dg{a) 


d In a 


^rn'^fRR{yH+g{a) + xa ^)Jd\na \3m'^fRR{yH+g{a) + xa 

3d+^Th^d2{l + 13) ^ 5 (a)(l - 3u;)^ ^(1 -b w)g{a) + fya"^^ -b {Jr - 1) -b 65 (a) - 6xa~'^ + 


din a 


fRRm^iyH + g{a) 4- ya 4) 


= 0.(71) 












































Taking into account the relations 
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Ti— “ 

a In a az 


( 72 ) 


din 


x2 d^ _ , d 

(1 + z) —^ + (1 + z) 
dz 


dz ’ 


we express the equation (71) in terms of the redshift, as follows 


(Pyn 

dz^ 


, j/i dyn ji/ jn „ 

+ —I" J2 yn + ^3 — 0, 


(73) 


(74) 


where 




1 


-3 + 


_ / fl -2 _ 1 

Qm'^fRRiyH + 5(z) + x(l + zf) \ ’ 


(75) 


T" — 

J 2 — 


{i + zy 


Ir 


'irn?fRR{yH + g{z)+x{l +z) ). 


(76) 


r// ^ d'^gi.z) _ 3 dg(z) _1_ 

^ dz2 (1 + 2 ) dz zfQrri^fRiiiyHg{z) + x{^ +z)^) 

X(1 + zf'^ + (77) 

/i{ ^ — (1 + z)-^^ + ‘2g{z) — 2x(l + z)^^ + + 2(;(z) + 2x(l + z)^ + 3^771^^ ^g(z)(l — 3 uj)^ 

Eq. (74) characterises the equation to be used for describing the cosmological evolution of the dark energy in the 
universe filled with collisional matter and radiation. We also solve numerically this equation for a given g{a) from 
(63) and we present the cosmological evolution of the dark energy scale yn = of collisional matter {w = 0.6) in 
comparison to the non-collisional matter as functions of the redshift. 


3^m^^2(l + 13) ( 3 (z)(l - 3u;) ) ( 5 (z)(l +w) + 



Figure 5: Cosmological evolutions of yn as functions of the redshift z for b = 0.19, Rs = 1.5 of collisional matter 

(blue) and the non-collisional matter (red). 
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In this rubric, we plot the evolution of the parameter of equation of state for dark energy wde{z) = Pde/pde, 
given by 


WDEiz) =-l + hl +— (78) 

3 vh dz 

derived by the continuity equation 

PDE + +'wde)pde = 0, (79) 

for collisional matter (w = 0.6) in comparison to the non-collisional one. 

As performed in the previous section, we are also interested to the evolution of others parameters. Theses parame¬ 
ters, the Hubble parameter H{z), the curvature scalar R{z) and the parameter Weff of the effective equation of state, 
can be described by still considering that the universe is filled by collisional matter and radiation. In the same way, 
we find the numerical solution for yn (74) for both collisional and non-collisional matters, and compare them. These 
parameters can be expressed as follows 


H{z) = \JTn?{yH + g{z) + x(I -h zf) 

obtained from equation (64), 


R = 3m^ ( 42/// -h 4g(z) - (1 -h - (1 + 


obtained by combining the equations (70) and (72), 


= -1 + 


2{l + z) dH{z) 
SH{z) dz 


(80) 


(81) 


(82) 



Figure 6: Cosmological evolutions of ujde as functions of the redshift z for b = 0.19, Rs = 1.5 of collisional matter 

for w=0.6 (blue) and the non-collisional matter (red). 
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Figure 7: Comparison of the Hubble parameter H{z) over z (left) and of the Ricci scalar R{z) (right), for h = 0.19, 
Rs = 1.5. The red line corresponds to non-collisional matter while the blue corresponds to collisional matter for 

w=0.6. 



Figure 8: Cosmological evolutions of Weff as functions of the redshift z for 6 = 0.19, Rs = 1.5 of collisional matter 

for w=0.6 (blue) and the non-collisional matter (red). 

For the Figs. 5, 6, 7 and 8, curves traducing the evolution of yniz), wde{z), H[z), R{z) and Weff{z), present some 
values strongly consistent with the cosmological observational data. These features prove that the fact of considering 
the coexistence of the collisional matter and radiation do not change the well know behaviour the different parameters, 
and conhrm that such consideration may be made. 


VI. CONCLUSION 

In this paper we undertake the f{R,T) theory of gravity. We focus our attention on the cosmological evolution of 
the parameter of deceleration q(z) and also the parameter Weff of the effective equation of state. The spacial aspect 
in this paper is that, besides the usual ordinary and the dark energy, new matter is assumed to participate to the 
matter contribution of the universe. This matter is self-interacting and called collisional matter, but with positive 
pressure. What we now generally is that the ordinary matter interacts with the dark energy, but any of them does 
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self interacts. From the cosmological considerations, its quite reasonable to introduce such a matter, i.e, a collisional, 
and search for its effects on the cosmological evolution of the universe. 

To this end, we considered two f{R,T) models as a sum of two gravitational models fi{R) and f 2 {T), respectively 
two functions depending the curvature scalar R and the trace T of the energy-momentum tensor. The general form 
of the f{R, T) models is constrained to the conservation of the energy leading to the form (32). The task has been to 
choose a suitable expression for the function f{R). Therefore, we considered, first, a generalized power-law in R, and 
the second expression is an exponential form in R. We plot the evolutions of q{z) and Weff for both the collisional 
matter, non-collisional matter and the ACDM and compare them. Our results show that, depending on the type of 
the model under consideration, the curve characteristic of the collisional matter approaches the ACDM one that the 
curve that describes the evolution of the dust. In the same way, the curve representative of Weff shows that its values 
are consistent with the observational, confirm the importance of considering the collisional matter in the study of the 
evolution of our universe. 

Moreover, we point out other type of matter, the so-called cardassian self-interacting matter, possessing a negative 
pressure and its corresponding parameter of equation of state has to be different from —1. This means that if such a 
matter is included in the content of the universe, there does not have vacuum energy. In this case, the same analysis 
have been done, as in the case of collisional matter, comparing the curves representative of this later to those of the 
dust and also the ACDM. Our results show that, for both the first and the second f{R,T) models, there is quite 
reasonable to consider the existence of the cardassian matter since the curve characteristic of this latter approached 
more the ACDM that the dust for some values of the redshift, confirming also that this kind of matter may be 
considered as a component of the universe, because also in this, the parameter Weff takes values strongly consistent 
with the observational data. 
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Appendix A: The energy-momentum tensor for a perfect fluid 


Let us assume that the Lagrangian density Cm is a sum of the kinetic term CmK and the one linked 
to the internal energy namely CmU^ i-e, 


Cm — CmK T CmU • 

Therefore, the energy-momentum tensor reads 


T — _ 


^ ("x/ "gCmK) ^ (x/ ~gCmu^ 








(Al) 


(A2) 


The well known expression of the kinetic energy-momentum tensor is expressed as T^v{k) = —pUfiUy. 
Let us consider CmU as 


i^mU = f{\Uk) = - 1^ . 


(A3) 


Using the method of Lagrange multipliers, the equations for the parameter A and the variables Uk are 
given by 


- 1 


and 


By setting = 0 and = 0, one gets ^g^'^u^Uv = 1 and X = p. 
Therefore, one can rewritten CmU as 


CmU = -P + p(^yJg^''U^Uy -1^ = -p 


(A4) 

(A5) 

(A6) 
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The internal energy-momentum tensor can be developed in the following from 

r7i _ f ‘^dCmU 

^ IJ.v{U) — 9iiv>-mU 


(A7) 


The second term can be expressed as 



(A8) 


The pressure p does not depend on the metric tensor , such that = 0, and then reduces 

to 


dCmU 1 


(A9) 



and finally, the internal energy-momentum tensor becomes 




(AlO) 


By using the contribution of the kinetic and the internal energy-momentum tensors and absorbing the 
sign (—) into the total stress tensor, one gets the energy-momentum tensor of matter fiuid as 


TfJ.1, = {p+p) u^Uy - pg^y. 


(All) 


It comes from the above calculus that the expression of the energy-momentum tensor is a consequence 
of the fact of considering that the pressure does not depends on the metric tensor g^''. 
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